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We study 2D quantum gravity on spherical topologies using the Regge calculus approach with the dl/l measure. 
Instead of a fixed non-regular triangulation which has been used before, we study for each system size four different 
random triangulations, which are obtained according to the standard Voronoi-Delaunay procedure. We compare 
both approaches quantitatively and show that the difference in the expectation value of R 2 between the fixed 
and the random triangulation depends on the lattice size and the surface area A. We also try again to measure 
the string susceptibility exponents through a finite-size scaling Ansatz in the expectation value of an added R 2 
interaction term in an approach where A is held fixed. The string susceptibility exponent 7g tr is shown to agree 
with theoretical predictions for the sphere, whereas the estimate for 7 str appears to be too negative. 



1. INTRODUCTION 

In the past few years Regge calculus has been 
extensively used in the study of quantum gravity 
PQ. In its usual form one studies regular simpli- 
cial triangulations of manifolds of a given topol- 
ogy, which are mostly hypertori. In this work we 
will investigate two-dimensional quantum gravity, 
where analytic calculations have shown that the 
internal fractal structure of the manifold depends 
very sensitively on the global topology. One uni- 
versal quantity is the string susceptibility expo- 
nent 7 s tr, which is the sub-dominant correction 
to the large area behavior of the partition func- 



tion Z{A) cx A^'- 3 e' 



-\ R A 



where Xr denotes the 



renormalized cosmological constant. The expo- 
nent 7s tr depends on the genus g of the surface 
through the KPZ formula 7 str = 2 - |(1 - g) gj. 

For the torus (g = 1) the Regge approach with 
the dl/l measure gives compatible results. For 
the sphere (g = 0) and topologies of higher gen- 
der, however, the situation is still unclear Pffij- 
One potential problem is that for the sphere only 
very small regular triangulations exist, such as 
the tetrahedron, the octahedron, and the icosahe- 
dron. In order to obtain triangulations on larger 
lattices, one either has to use non-regular trian- 
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gulations, with a few special vertices which might 
spoil the finite-size scaling (FSS) behavior, or to 
resort to random triangulations. In this work we 
use Monte Carlo (MC) simulations to study in de- 
tail random triangulations of a sphere and com- 
pare the results with our earlier results obtained 
by using the triangulated surface of a cube as 
spherical lattice 0. In addition we present es- 
timates of 7 s tr and the related exponent 7g tr |Hj 
on the random lattices, employing the novel FSS 
method of Ref. [H] ■ 

2. MODEL 

For a measurement of the string susceptibili- 
ties in Regge calculus one needs to introduce a 
curvature square term in the action, and then de- 
duces from its expectation value an estimate of 
7str and 7g tr through FSS analyses. We therefore 
considered the partition function 



Z(A) = Vfi(q)e 



-E, ^ 



-aRf 



6(£Ai-A), (1) 



where Rf = Sf/Ai denotes the local squared 
curvatures. The Ai are barycentric areas and 
Si = 2n — J2tDi are tne deficit angles, with 
0i(t) being the dihedral angle at vertex i. The dy- 
namical degrees of freedom are the squared link 
lengths, q = I 2 , which stand in a linear relation 
to the components of the metric tensor g. We 
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used the simple scale invariant "computer mea- 



sure" T>fi(q) = f| 



F e ({q i:j }). The nota- 



^fe) qij _ 

tion is identical to that used in Ref. 

The only dynamical term is the i? 2 -interaction, 
because we held A fixed during the update. The 
R 2 coupling constant a sets a length scale of ^fa, 
and A := A/ a can be used to distinguish between 
the cases of weak i? 2 -gravity (i> 1), where the 
KPZ scaling is recovered, and strong i? 2 -gravity 
(A < 1) where it was found |BJ, that Z(A) oc 
j^iU,-^ e -s c /A e -\ R A-bA^ witn tne classical action 

S c = 167r 2 (l — g) 2 , some constant b, and 7g tr = 
2-2(l-. 9 ). 

As global lattice topology we used a randomly 
triangulated sphere constructed according to the 
Voronoi-Delaunay procedure, see Fig.^for a sam- 
ple lattice. In this way we can control the influ- 
ence of non-regular triangulations. For spherical 
topologies we have the relations No — 2 = N 2 /2, 
N Q - 2 = iVi/3, and 2N X = 3N 2 , where N ,Ni, 
and N2 denote the number of sites, links and tri- 
angles, respectively. 



the necessary information to extract the string 
susceptibilities. A comparison with J2J yields 

d x (A) = b A + 7str - 2 + 0(1/ A) (4) 

for A ^> 1 and 

dx{A) = s c /A + j' stl .-2 + b Q A + o(A 2 ) (5) 

for A <c 1 . If we plot d\ versus A we expect to see 
a linear behavior for very large A, and a divergent 
behavior for small A, from which we can extract 
7str as well as 7g tr . Because R 2 in (PJ becomes 
infinite in the continuum limit N2 — ► 00, it was 
suggested to add a non-scale invariant part qfj to 
the measure and fine-tune a such that it cancels 
the divergent term do- However, a trial simulation 
showed within error bars no change in the relevant 
coefficient di{A) |7]. 




3. FINITE-SIZE SCALING 

All previously used methods to extract 7 str 0- 
E] are plagued by inconsistencies, as has been dis- 
cussed in our earlier work [7] . We suggested a new 
approach which is in spirit much closer to the con- 
tinuum analysis of Ref. [Hj as it uses A as the dis- 
tinguishing parameter between weak and strong 
R 2 gravity. The dimensionless expectation value 
R 2 := a(J2i R 2 ) can be shown to depend only on 
N2 and the dimensionless parameter A. Sending 
N 2 — > 00 one expects R 2 (A, 00) to be expandable 
in a power series, whose first three terms read as 




R 2 (A) 



+ b a A + b 1 +b 2 /A + 



(2) 



Figure 1. Randomly triangulated sphere with 
Nq = 1500. 



where b = —(\r - X)a, h = j stI — 2 for A > 1, 
and 6 = -b- (Xr - X)a, h = j' stI - 2, b 2 = S c 
for A < 1. 

Expanding R 2 (A, N 2 ) at constant A we obtain 

R 2 (A, N 2 ) = N 2 do{A)+d 1 {A)+d 2 (A)/N 2 + . . . (3) 

The next step is to expand the coefficients di as 
a power series in A. The coefficient d\ carries all 



4. SIMULATION 

For each lattice size we generated four different 
randomly triangulated spheres constructed ac- 
cording to the Voronoi-Delaunay procedure. Usu- 
ally, the size of the lattices varied from 218 up to 
17498 lattice sites, corresponding to 648 - 52 488 
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link degrees of freedom, or 432 - 34992 triangles. 
To update the links we used a standard multi-hit 
Metropolis algorithm with a hit rate ranging from 
1,. . . ,3. The area was kept fixed with a value of A 
in the range of 7 - 1 800. For each run on the four 
copies we recorded about 20 000 - 50 000 mea- 
surements of the curvature square R 2 := ^ R 2 
on every second to fourth MC sweep. The sta- 
tistical errors for each copy were computed us- 
ing standard jack-knife errors on the basis of 20 
blocks. The integrated autocorrelation time t R 2 
of R 2 was usually in the range of 5 - 10. As final 
error in the average of the four copies we used the 
standard root mean square deviation. 

5. RESULTS 

We first begin with a comparison of the raw 
data for R 2 /No obtained on the randomly trian- 
gulated sphere and our earlier data produced on 
the surface of a cube, see Fig. [21 The difference 
in R 2 depends, as could be expected, on the lat- 
tice size, such that the difference between the two 
triangulations decreases as Nq increases. This is 
also true for the difference in R 2 between different 
copies of the random triangulations. Noteworthy 
is, that for the larger system sizes the copies as- 
sume almost the same value within their statisti- 
cal error. 

However, the difference in R 2 also depends on 
A such that for small A the value of R 2 is larger 
on the cube, whereas for large A the value of 
R 2 is smaller on the cube than on the randomly 
triangulated sphere. This effect is not visible 
among the four copies of the random triangu- 
lations. Strictly speaking this shows that Z of 
eq. QJ and R 2 of eq. © depend also on the way 
the manifold is triangulated, i.e. the incidence 
matrix. Especially for large A, i.e. the region 
which determines 7 s t r , the values for d\ depend 
on the triangulation, and one can expect similar 
values only for very large lattices. 

We used linear two-parameter fits in eq. 
to extract from our raw data for R 2 (A, N2) the 
values di(A). A number of data points on the 
smaller lattices had to be discarded until the fit 
reaches a sufficiently high quality. All data points 
for g?i obtained in this way are shown in Fig. 
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Figure 2. Comparison of raw data for R 2 /No ob- 
tained for various lattice sizes on the non-regular 
cube and the randomly triangulated sphere. 



In the next step we fitted d\ (A) for small values 
of A < 120 according to the Ansatz (JSJ, which 
yields S c = 162(4) and 7^ tr = -0.2(4) with a 
total x 2 — 2.4, see Fig.0] This result is perfectly 
compatible with the theoretical prediction S c = 
16tt 2 « 158 and 7 ^ tr = 0. 

For large A we can employ Ansatz How- 
ever, there are only three data points with a suf- 
ficiently large A available. The linear fit yields 
7 str = —10(2) with x 2 — 0.8. It is not clear, how- 
ever, if we are already in the asymptotic regime, 
which might set in at much larger A. A further 
difficulty is that one is not interested in the slope, 
but in the intersection of the fit with the y— axis, 
which is far from the location of the points used 
in the fit. The systematic uncertainty on our esti- 
mate of 7 s tr is therefore hard to estimate, but it is 
interesting to note that our value is too negative, 
which is just opposite to what has been claimed 
in Ref. ^ by using a different FSS method. 

6. CONCLUSIONS 

Random triangulations appear to be a good al- 
ternative for topologies where no large regular tri- 
angulations exist. They show good scaling behav- 
ior, and the differences between different copies 
of the same area decrease as the system size in- 
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Figure 3. Results for di(A) vs. A together with Figure 4. Blow-up of Fig. 3 for A < 120 yielding 

fits in the regimes A < 120 and A > 400 yielding an estimate for j' str . 

estimates for 7^ tr and 7 s t r - 



creases. In this way they can provide a "typi- 
cal" lattice for the evaluation of expectation val- 
ues with the partition function of eq. @J. 

The quantitative difference in R 2 between the 
non-regular triangulation and the random trian- 
gulation of the sphere depends on both, A and 
No. In this way one will obtain on the usually 
used system sizes different values of di(A). The 
difference seems to be negligible for small values 
of A, so that 7g tr can be consistently obtained 
on both, the cube and the randomly triangulated 
sphere. However, the difference becomes impor- 
tant for large values of A, and is thus a potential 
problem for the determination of 7 s tr- 

Our FSS method of fitting at constant values 
of A gives results for 7g tr which are compatible 
with the theoretical prediction. In contrast to 
Ref. we employ a consistent FSS scheme and 
also much larger lattices. It would be interest- 
ing to test if contrary to [H] also for topologies of 
higher gender the theoretical expectations for 7g tr 
can be confirmed. 

Due to the fewer data points, only a crude es- 
timate for 7 s tr could be obtained which, however, 
appeared to be too negative compared to the KPZ 
theory. This is exactly opposite to what has been 
found in 0] with a different FSS Ansatz. We at- 



tribute this discrepancy to their method, which 
in our opinion [7j, bears conceptual problems for 
large values of A. It is yet unclear, if our system 
sizes are already in the asymptotic scaling regime, 
so that the potential danger of systematic errors 
is still very large. 
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